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A flexible theoretical model plasma which can be deformed to fit most measured elec- 
tron-ion-altitude profiles is employed together with available geophysical data on the iono- 
sphere to evaluate reflections and transmissions during quiescent and disturbed propagation 
conditions. The reflections and transmissions in the ionosphere are determined rigorously 
with the aid of the classical magneto-ionic theory. The complex indexes of refraction of the 
medium are deduced, and a coupling in the plasma between ordinary and extraordinary, up- 
going and downgoing modes of propagation is investigated. The corresponding reflection 
and transmission coefficients are then calculated, and certain phenomena which can be ex- 
pected as a result of the action of a solar disturbance on the reflection process are predicted. 

The disturbance of solar origin, investigated as an application of developed techniques, 
influences the reflected and transmitted LF waves in the lower ionosphere in a complicated 
manner. However, the high absorption phenomena exhibited by high frequencies do not 
seem to exist for the ])lasnia profiles investigated with the classical magneto-ionic theory. 

The electron-collision frequencies of the classical magnetic-ionic theory are modified 
to introduce a collision frequency proportional to the electron energy, and the changes 
necessary in the formulation of the classical theory as a result of such a consideration are 
presented. 

1. Introduction 

Low-frequency radio navigation and communication systems utilize or are capable of 
utilizing ionospheric modes of propagation, especially those systems with great distance trans- 
missions. Indeed, the nature of low-radiofrequency wave propagation around the earth is, 
iu large measure, determined by the shape of the lower ionosphere electron-ion density transi- 
tional region, i.e., the lower 40-100 km electron density-altitude, N{h), collision frequency- 
altitude, v{h)j profiles of the ionosphere. Previous theoretical treatment of the reflection 
process [Johler and Walters, 1960] at low radiofrequencies uses a sharply bounded model 
electron-ion plasma. 

Experimental evidence [Johler, 1961] indicates that such a profile is not always a valid 
model. Thus, especially in the case of an ionosphere disturbance, it seems to be appropriate 
to utilize a theory which takes into consideration in more detail the changes within the iono- 
sphere which affect the reflection process. This paper extends the theory previously pre- 
sented [Johler and Walters, 1960] to a flexible theoretical model which can be mutilated to fit 
rigorously almost any measured iV(A), v{h) profiles. The theory can therefore treat physical 
changes within the ionosphere as such changes affect the reflection and transmission processes 
and, in particular, the theory can treat a plasma reflection and transmission region with elec- 
tron-ion density transitions of quite arbitrary form. 

The newly developed model and associated techniques of analysis are then applied to 
certain geophysical and theoretical data in the form of the iV(A), v{h) profiles of the lower 
ionosphere during both quiescent and disturbed propagation conditions. The effect, if any, 
of a solar disturbance upon communication circuits which utilize low-frequency propagation is 
ascertained. 



1 This is an extended and revised version of a paper presented at a special meeting on disturbances of solar origin, of the Avionics Panel of the 
Advisory Group for Aeronautical Research and Development (AGARD), NATO, Naples, Italy, May 1961. 

2 This work was sponsored by the United States Air Force, Rome Air Development Center, GrifEss Air Force Base, N.Y,, and performed in 
connection with tasks 4 and 5 of delivery order AF 30(602)-2488 dated March 29, 1961. 
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The classical magneto-ionic theory is modified by the introduction of the concept of an 
electron-collision frequency proportional to the electron thermal energy [Jancel and Kahan, 
1955; Dingle, Arndt, and Roy, 1959; Molmud, 1959; and Phelps, 1960] into the analysis of the 
reflection and transmission process at the plasma with transitional electron-ion density, 

2. Determination of the Indexes of Refraction of the Medium 



The lower boundary of the model electron-ion plasma, figure 1, below which (^<C0) the 
ionization is nil (N=0), is considered to be the x^-plane. The region above the X7/-plane 
(0>O) is characterized by an electron density, N, and a collision frequency, v, which vary 
arbitrarily with respect to altitude, z. A plane wave, Ei, 



Ei= \Ei\ exp i(o:t---rjDj L 



(1) 



is incident upon the plasma from the lower region (2<C0), where the index of refraction (^<C0), 
'n=Vo='^) and c is the speed of light, c~3X10^ meters/second. The quantity D can be defined 
with respect to quantities in figure 1, 



D=x sin ct)i sin 0a+2/ sin 0^ cos 0^+^ cos (^j. 



(2) 



The field, Ei, is varying harmonically in time, ^, at a frequency /==co/27r. The field also is 
assumed to be incident upon the xy-ipil^me at an angle of incidence, 0^, and a direction of propaga- 
tion or magnetic aziniuth, 0^, measured clockwise from the ^^ -plane, figure 1. The earth's 
magnetic field vector, H^, is contained in the ?/0 -plane at a dip or inclination angle, / (measured 

from the horizontal). 

A resultant wave, Et, transmitted into the model ionosphere (^>0), figure 2, is then 
assimied to have the form 



where in the model plasma 



Et=\Et\ exp U {c^t-'^rjD\\, 



r]D=x sin 0i sin 0a+2/ sin </)j cos 0a+;^f, 



(3) 



(4) 





Figure 1. Coordinate systems. 



-ORDINARY (o)MODE 
_^ - EXTRAORDINARY ( e ) MODE 

' "Pi ! UPGOING (i) AND DOWNGOING (r) 

WAVES ARE COUPLED AT EACH BOUNDARY 

Figure 2. Structure of the flexible plasma model, 
illustrating ordinary and extraordinary modes of 
propagation, coupled at the boundaries. 

Each Zn, 71=1,2,3 ... p, becomes smaller as the number of slabs, 
p, is increased, until a stable reflection is obtained. Zp is always 
infinite, where each A^'n, vn corresponds to median values of the in- 
tervals, Zn, respectively of the particular profiles under investigation. 
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in which f is in general a complex number, the value of which will depend upon altitude, s, 

r=f(2). 

The quantit}^ f is determined by a simultaneous solution of Maxwell's equations,^ 

VX^+Mo^=0, (5) 

__ _ ?\fr 

VX^-J-eo^ = 0, (6) 

and the equation of motion of the electron, 



m'^+mgV+^oe (\'Xl^n?)+eE=0, (7) 



with the electric field, £',_volts/meter ; the magnetic field, H, ampere-turns/meter; the convection 
current density, J=—NeV; iV electrons per cubic meter (electrons/cnrMs frequently employed 
to describe A^ with charge, e; mass, 7n; and vector velocity, V, where ^mq and eo are the permea- 
bility and permittivity of space respectively. The constant real collision frequency, g^v, of 
the classical magneto-ionic theory is employed initially, but will be generalized to a complex, 
frequency-dependent parameter, ^, during the course of the paper. 

The vectors V and II can be eliminated, whereupon it can be concluded that the field E 
exists in the medium if a quartic m f is satisfied [Booker, 1934, Johler and Walters, I960]* 

ad'+ast'+(i2^'+ctd+ao=0, (8) 

where 

a.^(,u,,,-ir [l-^-^J+(sin^..-l) [1+^,+-^]+-^, (9) 

«^=2^P^(smV.-l). (10) 



^-{2t-CTj+^)}(-^^<^-^^- 



fiTCti I s—2 1 



(11) 



hJiraL ^ ^^^2 



a^=2 ^^^,_^,. -ai sec' 0^, (12) 

where v=Vn, the collision frequency of the particular slab, n, figure 2, under investigation,^ 

(15) 

(16) 

(17) 

ar=sin (^j sin </)a, (18) 

(19) 
(20) 



/l= 


COffCC 




K= 


=— A sin /, 




hr- 


=h cos /, 




aL^ 


=sin0< cos 0a, 


aT= 


a^ 


=sin (j), 






0^N^ = 


--Ne'/Km, 



3 Rationalized mks units are employed. 

* See eq (21) of paper by Johler and Walters [I960]. 

«In previous papers [Johler and Walters, 1960; Johler, 1961a] 00^ = ^//, ^cr=o)N' 
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where N=Nn the electron density of the particular slab, n, figure 2^ under mvestigation, 

o3H=iJLoeHni/m, the gyrofrequency,^ (21) 

K=l/cVo=6o. (22) 

Two pairs of roots, f , can be found by previously described methods [Johler and Walters 
1960] where each root represents either an ordinary or an extraordinary^, upgoing or downgoing 
wave propagated in the model plasma. 

The detailed structure of the flexible plasma model is illustrated, figure 2, as a stack of 
plasma slabs of arbitrary thickness (but consistent with computation efficiency) Zn except the 
topmost slab of thickness Zp= co . The number of such slabs, p, is also quite flexible, since the 
notion is implied in this analysis that the measured electron density-altitude and collision 
frequency-altitude profiles can be approximated to any desired accuracy by decreasing z^ and 
increasing p simultaneously until a stable reflection process is obtained. 

A constant electron density and collision frequency with respect to altitude, z, is of course 
assumed for each slab, z^; and associated with each such slab a set of four roots, ^=^n, is found 
to exist. Two of the roots will exhibit a negative miaginary, Im f negative, corresponding to an 
upgoing wave (+0 direction, figs. 1 and 2) and two of these roots will correspond to a positive 
imaginar}^, Im f positive, corresponding to a downgoing {—z direction figs. 1 and 2). Except 
for the topmost slab, it is necessary to consider both upgoing and downgoing waves in this analy- 
sis, whereas in previous analysis [Johler and Walters, 1960] only the upgoing waves were 
considered. The previous analysis implied that the upgoing waves were completely transmitted 
or absorbed before another reflection could occur, and all reflections occurred at the first bound- 
ary. Such a model must be applied with considerable caution. But the restrictions of such a 
model are removed by this analysis and, indeed, the sharply bounded model can be replaced in 
the limit by a continuum in which the electron density and collision frequency vary with alti- 
tude as a smooth but quite arbitrary function. This of course permits the use of the detailed 
measured profiles of the lower ionosphere. 

The treatment of a continuously varying medium by approximating the medium with one 
or more slabs of uniform composition has been utilized by many authors [Hines, 1951 ; Ferraro 
and Gibbons, 1959; Brekhovskikh, 1960; Wait, 1960 and 1961]. Such methods are exploited 
in this paper and carried to the limit, such that the number of slabs, p, for each calculation 
depends upon the computation precision required and the particular values of the electric and 
geometric parameters. 

It is, however, necessary to distinguish between ordinary and extraordinary modes of 
propagation for both the upgoing and downgoing waves. This is accomplished by an examina- 
tion of the form of the index of refraction function with respect to frequency and altitude (or 
electron density and collision frequency which varies with respect to altitude). Thus, the index 
of refraction, rj (as defined by eq (4)), 

rj'=t'+sm'^„ (23) 

is detailed for each frequency and slab Zn, ^1=^71, and the upgoing ordinary and extraordinary 
^1^0.6 and the downgoing 7]^^^.g ordinary and extradordinary function continuity is examined 
in detail as a function of frequency to determine the crossover point of the functions^ for 
each slab or electron density under consideration. The Re rj was employed in this analysis. 
Thus, above the crossover point the ordinary wave was considered to be the greater of the two 
roots Re T/^^Re Ve and below the crossover point the ordinary wave was considered to be 
the lesser of the two roots. Re r/o<Re rje. At precisely the crossover point, the two, of course, 
are identical, r}o='ne which point, for examples with low electron density examined in this 
analysis, was found to be below the critical frequency, /at =co^/27r. The absolute distinction 
between the two modes of propagation remains quite arbitrary as in previous analysis [Johler 
and Walters, 1960], but the analysis must be consistent between each slab and consistent 
within each slab for upgoing and downgoing waves. 

6 Index of refraction, rj, continuity analysis to determine crossover was coded as a logic routine for the CDC- 1604 computer .as was the whole- 
of this analysis. 
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the E field. 

*'x' 1 *^°' *a * ''y' 1 '^°' *i *'" *a * ^x' r '^°* *a ' ^y' r '^°^ *i ^'" *a ^ '^io ^yio * ^"'le ^yie * ^ro ^yro * ^re *'yre 

q|1> e('>^ exp I -i ^ ,j t|^> 1 . Qf> E<'.>^ exp ( -. ^ .^ ^|^» ) . Q^^^ e''^^ exp [ -, ^ z^ t,['^ ] . Q^'j E^^*^ exp [ -. ^ z, ;<'» ] = q\'J E<^^ . q|^> eJ^> . o'^j e'^^^ . q['J ^[\[ 

q|^' E<^> exp [ -i ^ z^;<^>) > Q<^> E<^> exp [ -> ^ z^ r.'fj ] , Q<^j E^^^^ exp [ -. ^ z^ ^['^ ] . q'^^E*^;^ exp [ - ^ z^ ;<^J ] = qJ^^* E^^> . q|^^' e'^^ . q['J z[\[ . o'^j e'^^^ . 

(24) 

Q<P-^' E<P-^' exp [ -i ^ z , t<P-" 1 i Q<P-'> E<P-" exp [ -. ^ z , C^P"^*) + q'P"^' e'P'^* exp ( -. i^ z , ;(P-1>1 . Q<P-'> E<P-1' exp [ -. ^ z , ;<P-'>1 = Q<P' E<P' . q'p' E<P> 
lo yio "^ ' c p-1 *io ' le yie "^ ' c p-1 le ' ro yro '^ ' c p-l ro ' re yre r^ i ^ p.i »re ' lO yio le yie 

the Ey field. 

r „ * * F A * F * F .„*.„.* F<'> * r<l> * F<'> + F<>> " 

x' i *a y- 1 *i *a " x' r *a " y' r *i "^°*«*a " yio yie yro yre ' 

yio P ' r 1 ^lo ' yie P ' F 1 ^je ' yru P ' c 1 ^ro ' yre P ' F I ^re ' " yio yie yro yre ' 

E<^> exD ( -i - z t'^'l + E*^* exQ I -i !^ z I.<^'| t E<^' exo [ -i !^ z t^^'l + E<^> exo t -i ^ z t'^'l = e'^' + e'^' + E<^' + E<^' 
yio *^ ^ c 2. io ' yie " '^ c 2 ie ' yro " '■ c 2 ro ' yre " ^ c 2 re ' yio yie yro yre 

E<P;'> exp [ -1 ^ Zp_, ;|P-"l . E«P;'> exp [ -. ^ Zp_, ;|P-'»1 . E<P-j> exp 1 -i ^ z^., ;^P-^'l . E<P;" exp [ -, ^ z^., C^P"") = E<P> . E^P> ^ ' 

the H field, 

E^. . cos ♦. sin ♦^ - E^, . CO. ♦^ - E^. ^ cos 0. sin ♦^ - E^. ^ cos ^^-[.^ p|J,> - ;<),>1 e(\>^ > [ a^ P<^» - t|^>] e'^, ^ ( a^ P<^j - ;<'j] e''^^ . (a^^ P^'j - c'/j] E<t>^ . 

l»L PJo* - ^!o'l <U «P I -' F ^1 '•lo'l ^ [-L ^le' - ^!e*l ^'.e "P I - F ^, ^!e'l ' l*L ^to' " ^I'o'l ^J'ro "P 1 - ? M ^o 1 ' 1»L ^te' ' ^rV' ^fre "P ' - F ^ C'l = i*L ^lo' " ^!o'l ^'1 ' [»L ^If " ^If 1 ^'i ' 1*1. ^it' ' 4o 1 ^r'o ' l*L ^le ' ^I'l ^r'e ' 

[*L ^!o' - ^!o*) K'i'o "P [ -' F ^2 4'*1 * t^'L PJf' - ^!e*l ^^li ^P I - F ^2 ^le'l ^ ^'l ^ro " ^[t'l ^ro "P 1 - F ^Z ^rl'l ^ ''l ^I'J " ^tt'l ^re "P I "' F ^2 ^'re'l^ l*L p!o' " ^If') ^li ^ 1*1. ^le' " '■[e^ ^[l ' t^U '^to " 4o ) ^r'o ^ t*L ^re " ^re'l ^^r'e ' 

00 

01 (26) 

l*L PJo"^' - ^*r**l *^l^;" ^"P [ -' ? -p-l ?.!r''j ^ I-L ^Ir" - ^!r"l =1^;" «? l - F ^.l ^!r"l ^ 1»L ^to"" - ^to'") ^t^/o' "P ^ - F ^ .1 ^iV^^ * [»L ^'re"" - ^'e'") ^^'^/e" ^P l - F ^ -1 ^re""l =1*L ^Io' - ^!o'l *^*^o ' t^L^le* ' ^le'l ^^"^1 

Uie H field. 
1 

E^. ^ cos ♦, cos ♦^ . E^. , sin .^^ - E^, ^ cos «. cos ♦, + E^. ^ sin ♦^ = [ ;<^' Q<^' - a^ p|''j eJ'^* . [ ;<^> Q<^' - a^ P<^>] E^',^ > [ ;^'j Q^'j - a^ P^'j] e''^^ . [ t,[[^ q'^ - .^ pf'J] e''^ 

1 ^o ""[o - *T ^!o*) ^t'lo "P 1 -' F ^ ^!o*J * 1 ^!e* '^le' " *T "le'l ^\\l "P ' -' F ^I ^le'l * ' ^to '^ro " »T ^lo > ^i'r'o "P ^ ^ F ^1 ^ro 1 ' ^ ^'re' "l'.* " 't ^'rl'l ^fr'e "P ' "' F 'l ^e > = I ^If' '^lo' ' *T ^lo'l ^li * I ^e' ^if ' *T ^l:'' ^i'i ^ I ^ro '^ro " *T ^^to' ^rl * ' C ^re ' »T ^re'l ^I're 
I ^*o' Q!o' - ''T p!o'1 =y1o -P 1 -^ F ^2 ^!o'l ^ 1 ^'J "le' " *T ^!e') ^^li "P ' - F ^2 4") ' i ^o '^ro " »T ^to > ^r'o "P I "' F ^2 ^ro 1 ^ I C ^re " >T ^'e 1 ^tr'e «P I "' F ^2 ^It'^ I ^''^lo' " ^T ^lo') ^\]i * ^ ^!e' ^le' " *T ^le'^ ^y!i ' I ^o '^'o " *T '^lo i ^i'ro ^ I 4t' "rf ' *T ^te'l ^'rl 

[ ,^(P-1) q(p-1) . ,^ p(p-l), ^(P^n exp [ -1 ^ zpj ;<P-1)) . [ ;<P-" Q<P-'> - a., P<P-1>] e(p^'' exp [ -. ^ z^^j ;|P-"] + [ ;<P;^' Q^P-" - ..^ P^P-l'l E<P/J> exp ( -i ^ z^., ^IP-^'] . ( ;<?"'> ulP-" - a^ P;P-'>1 eJP;'* exp [ -, ^ x^^ ^P"'*] = tJP* q|p' - a^ pJP'l E^P> i [ ;|?> q|p> - a^, P<P»] £<?> . . (27) 
the H^ field. 

(28) 



3. Determination of the Reflection and Transmission Coefficients 

The reflection and transmission coefficients are determined b3' _tlie boundary conditions 
which express the principle of continuity of the tangential E and H fields and the normal 
H fields at each boundary, figures 1, 2, of the m.odel plasma. These fields are equated im- 
mediatel}^ above and immediately below each boundary. (See eqs. 24 to 28 on p. 85.) 
where 

Q=^> (29) 

P=§' (30) 

and the subscripts io, ie, ro, re refer to upgoing ordinary, upgoing extraordinary, downgoing 
ordinary, and downgoing extraordinary modes of propagation. The superscripts n—lj 2, 
3 . . . , p—l, and p refer to the particular medium, figure 2, under consideration, with cor- 
responding values Nj v, which together with the magnetic parameters, Hm, (t>a, I, and the angle 
of incidence, </)^, determine the index of refraction, (23). The process is of course continued 
until convergence is obtained on the particular profiles of electron density and collision fre- 
quency under investigation. Notice that the topmost slab, Zp=co , ^=i>j has only two terms 
representing the two upgoing waves. The downgoing waves do not exist since the topmost 
slab is assumed to be of infinite, homogeneous extent. However, this analysis can readily 
determine the transmission coefficient of a finite slab of plasma with an N{h)=N{z), and v{K) = 
v{z) profile which varies with altitude, s, simply by replacing the last two terms of eqs (24 to 28) 
with 

E^'t cos ^a+Ey^t COS (t)i sin 0«, 

—Ex't sin ^a+Ey^t cos 0^ cos 0„, 

Ex't cos </)i sin (j)a—Ey^t cos </>«, (31) 

Ex't cos (j)i cos (l)a+Ey^t sin (ji)^, 
— E';,/^ sine/),, 

respectively, where Ex^ty Eyt, are the components of the transmitted fields. 

The refiection, T, and transmission .T", U, coefficients can now be defined 

rp __ L^y'r rjif __ L^v't j-r(n) -^ yio jj (n) -^ yro 

-^ ee TTi -^ ee rp ^ eio jp ^ ero — 777 

^V'i ^v'i J^y'i JOjy'i 

rp __ L^x'r rpr __ J^x't ■rr(n) ^ yiQ Tj {n) - ^ yro 

•^ em Tp -^ em JP ^ mio jp ^ mro — tti 

J^y'i ^y'i -^ x' i J^ x' i 

rp _ J^y'r rp, __ J^y't Tjin) -^ yie rj (n) ^ yre 

-^ me JP -^ me JP ^ eie jp ^ ere — ^ 

■'^z'i J^xH J^y'i -t^y'i 

E , E .s TT (n) JP in) 

rp ^ J^x'r rp, _ J^x't jji^ri) _ J^ yie jj („) _ ^ yre /qqN 

■^ mm JP -^ mm jp ^ mie tp ^ mre — ttt W-"/ 

■^x'l J^x'i J^x'i -tl'x'i 

where, figure 2, n=l. 2, 3 . . . , y—l,/D, and the four T's are defined as a transmission coef- 
ficient into a topmost slab of infinite extent, s, and zero electron density, iV=0, rjp=7jo=l. 

The reflection coefficients. Tee, Tem, Tme, Tmm, describe the reflection into the region below 
the model plasma, 2<0. Thus, Tee refers to vertical electric polarization of the incident wave 
and a corresponding vertical electric polarization of the reflected wave. Tem refers to the 
generation of the abnormal component (horizontal electric polarization) by the incident 
vertical electric wave. Similarly, Tmm refers to vertical magnetic incident wave and vertical 
magnetic reflected wave, and Tme refers to a corresponding abnormal component (horizontal 
magnetic). 

The transmission coefficients, U, refer to transmission at the particular point in the iono- 
sphere under investigation to which the wave has penetrated. There are thus four modes of 
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propa|>:atioii and two polarizations and hence eight types of transmission coefficients: upgoing 
and (lowngoing, ordinary and cxtraorchnary, vertical electric and vertical magnetic. 

Tluis, the nature of the wave reflected from the plasma, the nature of the wave trans- 
mitted through the plasma, and indeed the nature of the wave progressing within the plasma 
are completely described by the coefficients determined by this analysis. 

Four of the boundary conditions (24 to 27) result in the matrix equation 

M, 1 *1. 2 *1. 3 *1,4 ^1, 5 ^1, 6 

*2, 3 ^2. 4 ^2. 5 *2. 6 *2. 7 *2, 8 *2. 9 *2, 10 

*3, 7 *3. 8*3.9 *3. 10*3. 11 *3. 12 *3, 13*3. 14 



P.P + 4- 



P.P + 9 



3, 11 3, 12 3. 13 "3. 14 



p. p + 4 • 



p. p + 9 



'^l. 1 ^1. 2 *^1, 3 ^^1.4 *^1. 5 '^l. 6 

^l, 3 *^2.4 ^2, 5 '^Z, 6 ^2. 7 ^l, 8 ^2. 9 ^Z, 10 

*^3. 7 *^3. 8 *^3, 9 ^3, 10 '^ 3. 11 ^3, 12 *^3, 13 ^^3. 14 



P.P+ 4- 



P.P + 9 



p. p + 4 • 



p. P + 9 



ee' 


^n.e 


(1) 


,,(1) 


eio' 


mio 


(1) 


n<^> 


eie' 


mie 


IVo 


yd) 

mro 


(1) 


,/l) 


ere 


mre 


(2) 


n(2) 


eio' 


mio 


(2) 


,,(2) 


eie' 


mie 


(2) 


„(2) 


ero 




(2) 


„(2) 


ere 


mre 


(3) 


h(3) 



u{p-l)y(p-l) 

eio' mio 

^j(p-l)y(p-l) 

eie' mie 

u(p-iVp-^) 

ero ere 

^j(p-l)^j(p-l) 

ere" mre 

u<P> u<P) 

eio' mio 

u(p) u(p> 

eie mie 



o, e \ o, m 



(33) 
Rearrangement of the matrix, such that zeros do not appear on the diagonal, ai,i to 
r/p,p+9, results in a matrix which can be solved by the previously described methods [Johler 
and Walters, I960]. (Note that the second subscript of the last term of each group is ;p + 9 
instead of f+ll since the topmost slab contains only the ordinary and extraordinary upgoing 
waves.) 



*1. 1 *1.2*1, 3*1.4 *1. 5*1.6 
''l.l '*1.2**1.3^1.4^1. S'^l.fe 
•^1. 1 *^1.2 '^l. 3 '^1.4 *^1. 5 *^1. 6 
*^1.1 **1, 2^*1. 3*^1. 4*^1. 5*^1.6 

*2, 3 *2. 4 *2. 5 *2. 6 *2. 7 *2. 8 *2. 9 *2. 10 
^Z. 3 ^1. 4 ''2. 5 ^l. 6 ''2. 7 **2, 8 ^Z. 9 *'2. 10 
''2. 3 ""2, 4 '^2. 5 *'2, 6 *^2. 7 ^Z. 8 *'2. 9 ^"2. 10 
**2. 3 '^2. 4 '*2. 5 **2. 6 **2. 7 *^2. 8 **2. 9 ^^2. 10 

*3. 7 *3. 8*3. 9 *3. 10 *3. 11 *3. 12*3. 13 *3. 14 
^l.l^l. 8 H 9 ''3. 10 ^3, 11 **3, 12 ^X 13 ^\ 14 
*^3,7 ^3.8 *^3. 9 *^3, 10 ^^3. 11 '^ X 12 "^S. 13 "z, 14 
'*3.7*^3.8**3.9**3. 10 ^3.11 **3. 12 **3. 13 **3. 14 



p. P + 4 

V P + 4 

S.P + 4 

S. p . 4 



*p.p+ 9 

''p. P + 9 

S. P + 9 

%. P + 9 







em' mm 




*oe' *om 






^ee' ^me 




^e- ^m 






u(i) u(i) 

eio* mio 










u(i) u<»> 

eie* mie 




<*oe' -^om 






ero' mro 










ere' mre 
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(34) 
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(p)i 



aoe = COS (l)i sin (l)a 
boe = COScl)i COS (f)a 
Coe^ — COS (l>a 
doe = Sm(l)a 



aom=<^OS <j)a 

Co7/i=cos (pi sin 0a 

^om=COS0t COS</)a. 



The analytic expressions for the complex numbers, Pl2^ Fif , P^V, P^,?, Q^^?, Qlf , Ql'S ^ 
QS-e^ can derived from the definitions (29, 30) by a simultaneous solution of Maxwell's eq (5, 6) 
and the equations of motion of the electron (7) with the following result: 
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(34) 



(35) 



where the particular slab, figure 2, n=l, 2, 3 . . . p-1, p, under consideration is designated 
by the notation, f =f^^\ P=P^'^\ Q=Q'^'^\ and, io, ie, ro, re refer to the four roots of the quartic 
in f for upgoing, downgoing, ordinary, and extraordinary waves respectively in the particular 
part 2n of the electron-ion plasma under consideration. 

The ratio of the fields Q^E^lEy and P=EJEy was introduced in the boundary conditions, 
eqs (24 to 28"^ as each slab was matched electrically to its adjacent slab such that the tangential 
E and H fields and the normal H fields were continuous across the boundary. Notice that onl}^ 
the first four sets of boundary conditions, eqs (24 to 27), were employed in the matrix, eqs 
(33 to 34) ; the fifth boundary condition served as a check on the entire computation, since it is 
automatically satisfied by the other four conditions. 



4. Monoenergetic Electron Collision Frequencies 

The recent work of Jancel and Kahan [1955], Dingle, Arndt, and Roy [1956], Molmud 
[1958], and Phelps [1960] has treated the collision process in an electron-ion plasma with greater 
rigor. It can be shown, as a result of these investigations, that the average collision frequency, 
g=v, in the equation of motion of an electron (7) is an effective value independent of the 
energy, u=^kT/e, where k is Boltzmann's constant, Tis the temperature, degrees Kelvin. How- 
ever, the investigations of Phelps [1960] indicate a strong dependence of the collision frequency 
upon the energy, u. 

Phelps has developed integrals of the electron velocity or conductivity tensor. The 
integration over a Maxwellian energy distribution, fo=[e/TkT]^^^ exp [—eulkT]^ involving a 
momentum transfer collision with gas molecules proportional to the energy u, v=v(u) is accom- 
plished for an electromagnetic wave of angular frequency, o:==2tJ. The magnetic field is 
assumed to be vertical. The results of Phelps are generalized in this paper to a magnetic field 
of arbitrary direction and introduced into the analysis of the plasma model with a transitional 
electron-ion density described in this paper. 

Indeed, Molmud [1959], demonstrated that the monoenergetic electron collision-frequency 
concept can be introduced into the equation of motion of the electron, 
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m^+mgV+nae(VXHJ+eE=0, 



(36) 



as a new collision parameter, g, which has replaced the classical; v. Thus, the classical theory 
implies, g^^v. 

The velocity, V, tensor can be written, 
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The coefficients of the tensor matrix can be reduced to the form 

(4) , (-1) 



Ar- 



A, 



A,= 



Br- 






(— co^eo), 



(— oj^^eo), 



(—ai%eo), 



(— co^-eo), 



_r (-cos^/) (-fsin^/) (-Isin^J) -] , 
D _r (sin /cos 7) (— I sin J cos /) (— | sin /cos /)"1 , ^ \ 

^ _r (; COS /) ( — 1 COS 1) 1 ,_ 2 , 

'lg{^0+i^r'g{ii2)+iii2A ^ '"'"''' 

^ _[ ~(sin /cos /) (— ^ sin /cos /) (— | sin /cos /) ~| ,_ , -, 
'~L y(no)+i"o + fif(0:)+iO: + r7(02)+'ifi2 J ^-"•^•*°^' 

(-sin^ /) (-1 cos^ /) , (-i cos^ /)' 



a 



{: 



,^(l]o)+^l^o ^(12i)+i^l2i ^(122)+^l]2 



] 
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(37) 



where 5i=— ^2, Ci=—Az, C2=Bzy and where 12o=co, 12i = co+co/y, 122=co— co^. A quantity, F, is 
found by integrating over the velocity distribution for monoenergetic electron-ion coUisions, 
v{u)^ at each frequency, 12=12o, ^1, ^2, 
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— u: 



[v{u)+ia\ du 



lirkTi 



3/2 



exp [—eu/kT]du, 



(38) 



which can be written in terms of Dingle, Arndt, and Ko}^ [1956] tabulated functions, ^5/2(^)9 
^3/2 W? 

(39) 



^=^ G 6) ^''^ \V~' \iv ^"' ©J ' 



whereT = v('u)/12,u=/:T/6, and j^('Z/) = akT/e=an. The complex collision parameter, ^, for the 
equation of motion can then be evaluated 



i=T-'-in= 



BeT-i[[mT + n\T\ 



(40) 



The quantity, g, is thus complex and there are three values corresponding to frequencies 
12 = Qo, ^u ^2. Furthermore, the value of g is also frequency dependent. The implication of 
the classical magneto-ionic theory that such a quantity can be represented by a real effective 
value without a frequency dependence is therefore not in general a valid concept. 

The complex collision frequency, g, can be introduced into the anah^sis by evaluating the 
coefficients of the quartic (8 to 13) as follows: 

ai=^—o)hl(C3+io:€o), (41) 

a,= -2c^hlaLBs, (42) 

a2=o:hl[A + B2+2Cs-al{B,+Cs)-aUA+Cs)]+icoeo[-A,a-B2Cs-AI+BI^^ 

(43) 

ai = -2ic^eoaL(A2As+AiBs) -2o:^e^aLB^{a:'— 1) (44) 

ao=AB2Cs-2A,A,B,+AIB2-ABI+AlCs+io^eo(A,C,+B2Cs+A,B2+Al+AI-Bl) 

+ ic0 6oaM-^lC^3-A52-^i-^D+^^6o4(-^2a-A^2 + ^i-^D+^'eg(-^l-52-Q 

+ cc'elaU2A + B2 + Cd+o^hl2l{A + 2B,+C,)+c^hlaWL(-A,-B2) 

+ c^'elaU-A,)-o:hlaiB2-ioohl(l-ay. (45) 

The functions fp (x) {p=7^^ k) can be evaluated by an application of quadrature techniques to 
the integrals [Johler and Walters, 1960], or other methods described by Dingle, and Roy [1956], 

h(x)=Q^}jj\Hx^+pT' exp (-p)dp, (46) 



Ji 



^'^=(¥)£ 



pi{x^+p'') 'exp {-p)dp. 



(47) 



The quantities P and Q are rewritten,^ with the aid of eq (37), by substituting for s in eq 
(34, 35) the corresponding values of g, %, 12i, 122, 






^coeoJ 
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l-a^ 



l-a2+^ 



C, 
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tweoJ L ^weoJ L ^'<'eoJ 



(48) 



(49) 



7 The constant eo or the factor im can be eliminated from eq (37, 41-45, 48, and 49) for convenience. The factor was retained so that eq (37), 
would be correct dimensionally, i.e., J= current, amperes per square meter. 
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5. Computations and Discussion 

The reflection coefficients, figures 5-24, were evaluated for the continuous electron-ion 
density transition, figures 3 and 4, with the aid of the classical magneto-ionic theory. The 
typical daytime- noon profile, figure 3 [Waynick, 1957], was selected to represent a quiescent 
ionosphere. The disturbed profile, figure 3 [Seddon and Jackson, 1958], was also selected to 
give an indication of the intensit}^ of the ionization during disturbed conditions. The region 
of the ionosphere below 1,000 electrons/cm^ was represented with a Gaussian el(H*tron density- 
altitude distribution [Budden, 1955], 

A^^A^max exp [-{z-2m^,fl¥l 

wliere the constants A^max, ^'max, ^md k/ were determined by the measured profile, > 1,000 
electrons/cnr\ In particular, .Vn,ax=26,000 electrons/cm.^ 2max--85,000 m, k' = 1 .92 X 10^ This 
profile was measured during a high-frequency ionosonde blackout observed at Ft. Churchill, 
Canada, 1332 CST, 15 November 1956. Data was therefore made available in the theoretical 
form of a continuous electron-ion density transition for a model plasma whicli represents 
approximate conditions of the ionosphere during disturbed phenomena. This ionosonde 
blackout is classified as a lo(^al auroral-zone blackout which is confined to the auroral zone or 
some segment thereof. 

The Nicolet/3 collision frequency [Nicolet, 1958; Sedden and Jackson, 1958] was employed 
in this analysis. Other profiles, vQi), are illustrated, figure 3 [Compton, et al., 1953; Gardner 
and Pawsey, 1953; Fejer, 1955]. These collision frequencies are applicable to the classical 
magneto-ionic theory employing the approximation, g^v in the equation of motion (7) of the 
electron. Recently Kane [1959 and 1960] re-evaluated the pro file; v{h), on the basis of mono- 
energetic electron collisions, v=v{u), from which the complex parameter, g, can be deduced. 
Furthermore, the reflection coefficients presented in this paper can be re- evaluated with the 
aid of such data and the theory of the monoenergetic electron collision frequency, which task 
will l)(^ I'eported separately for comparison with the classical theory. 

The computation of the field strength of a radio wave propagated about the earth is 
dependent upon the reflection process at the ionosphere [Johlcr, 1961 b and c] and the mode 
theory [Wait, 1960] or a geometric-optical theory requires the introduction of reflection co- 
efficients into the analysis. Thus, the essential nature of ionospheric propagation is described 
by the reflection coefficients, Tee, Tem, for vertical polarization, and Tmm, T^e for horizontal. 

The reflection coefficients, Tee, Tem, Tme, Tmm, were evaluated, figures 5 through 11. For 
the typical daytime-noon model ionosphere as a function of frequency, /=10 kc/s to 1,000 
kc/s, with various values of magnetic azimuth, (t>a=0, 45, 90, 135, 180, 225, 270, 315 deg; at a 
temperate magnetic latitude, 7=60 deg, assuming a typical magnetic field intensity; Hr,t=0.5 
gauss and grazing incidence transmission, (^^=81.79 deg (d = 1,000 statute miles for a single-hop 
geometric-optical transmission) . 

The reflection coefficients, | Ted, \Tmm\ show a steady decrease in amplitude as the frequency 
is increased. A small, perceptible rise is noted near the gyrofrequency, however. Such a 
steady decrement in amplitude could be imitated by the sharply bounded model ionosphere, 
[Johler, 1961 b and c] by suitable selection of a fictitious low-electron density. But such a 
model is a primitive theoretical representation of the true reflection process. 

The magnitude of the abnormal components, \Tem\, \Tme\ also show a decrement as a 
function of frequency. But these components are quite small (<0.1 at 10 kc/s, for example). 
The corresponding disturbed-blackout reflection coefficients are also illustrated. The angle 
of incidence, 4>i, was changed slightly as a result of a lowering of the ionosphere. The principal 
cause of change as a result of this profile is nevertheless the re-distribution of the N{h) profile 
from the quiescent to the disturbed conditions. It is interesting to note an increased atten- 
uation at the higher frequencies, 100 to 1,000 kc/s. However, the high attenuation or ionosonde 
blackout which characterizes high frequencies does not seem to exist. Indeed, an enhancement 
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Figure 3. Electron density-altitude 'profiles N(h) 
{quiescent and disturbed-blackout) employed in the 
analysis. 
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PiGUEE 5. Reflection coefficients {amplitude and 
phase) of the lower ionosphere for low frequencies ^ 
illustrating the frequency dependence of the classical 
theory (f>a=0, 180. 
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Figure 6. Reflection coefficients {amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the frequency dependence of the classical 
theory, <i>a=90, 270, vertical polarization, Tee, Tem. 
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Figure 7. Reflection coefficients {amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the frequency dependence of the classical 
theory, (f)a = 90, 270, horizontal polarization, Tmm, 



Figure 4. Collision frequency-altitude profiles v{h) 
for classical magneto-ionic theory. 

Nicolet/3 was employed in the analysis. 

in the field is anticipated from this profile at VLF «30 kc/s). Whereas there is no experi- 
mental evidence to establish the non-existence of ionosonde blackout phenomena at LF/VLF, 
a search of the literature did not reveal any report of very high attenuation of LF/VLF during 
such D-region events. 
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Figure 8. Reflection coefficients {amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the frequency dependence of the classical 
theory f 4)a=4^y 135, 
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Figure 9. Reflection coefficients {amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the frequency dependence of the classical 
theory, 4>a=225, 815. 
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Figure 10. Reflection coefficients {amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the action of the Lorentzian force {earth^s 
magnetic field) f=10 kc/s. 
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Figure 11. Reflection coeffiicients {amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the action of the Lorentzian force {earth's 
magnetic field) in the classical theory, f=20 kc/s. 
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Figure 12. Reflection coefficients {amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the action of the Lorentzian force {earth's 
magnetic field) in the classical theory, f= 50 kc/s. 
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Figure 13. Reflection coejpxients {amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the action of the Lorentzian force (earth^s 
magnetic field) in the classical theory, f=100 kcjs. 
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Figure 14. Reflection coefficients (amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the action of the Lorentzian force (earth's 
magnetic filed) in the classical theory, f=200 kc/s. 
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Figure 15. Reflection coefficients (amplitude and 
phase) of the lower ionosphere for low frequencies 
during blackout conditions illustrating the frequency 
dependence of the classical magneto-ionic theory. 

\ Tee\ and \Tmm\ are also illustrated for comparison with quiescent 
conditions (see fig. 5-9 for complete quiescent reflection coefficients) . 
00=0. 180. 
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Figure 16. Reflection coefficients (amplitude and 
phase) of the lower ionosphere for low frequencies 
during blackout conditions illustrating the fre- 
quency dependence of the classical magneto-ionic 
theory. 

I Tee\ and | Tmm\ are also illustrated for comparison with quiescent 
conditions (see fig. 5-9 for complete quiescent reflection coefl3cients). 
0a=9O. 270, vertical polarization, Tee, Tem. 



The effect of the Lorentz force, ixQe{Vy^Hm), (7) or the action of the earth's magnetic 
field, on the reflection process is also illustrated, figures 12 through 15. The reflection coeffi- 
cients are illustrated as a function of the magnetic azimuth, 0a. It is interesting to note that, 
as in the case of the sharply bounded model [Johler, 1961a], the propagation into the west, 
(/)a=270°, shows a smaller reflection coefficient than the propagation into the east, (^^=90°. 
Thus, there is clearly a non-reciprocity (interchange of transmitter and receiver) in the propa- 
gation which enhances the field propagated into the east relative to the field propagated into 
the west. However, the characteristic variation with respect to magnetic azimuth is consider- 
ably masked at the higher frequencies (>200kc/s) for the highly absorbing daytime-noon 
model ionosphere employed. 
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FiGUKE 17. Reflection coefficients {amplitade and 
phase) of the lower ionosphere for low frequencies 
during blackout conditions illustrating the fre- 
quency dependence of the classical magneto-ionic 
theory. 

\Tee\ and \Tmm\ afo also illustrated for comparison with quioscont 
conditions (see fig. 5-9 for complete quiescent reflection coefficients). 
0a=QO, 270, liorizontal polarization, Tmm, Tme. 
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Figure 18. Reflection coefficients (aniplilude and 
phase) of the lower ionosphere for low frequencies 
during blackout conditions illustrating the fre- 
quency dependence of the classical magneto-ionic 
theory. 

\Tee\ and \Tmm\ are also illustrated for comparison with quiescent 
conditions (see fig. 5-9 for complete quiescent reflection coefficients). 
<},a=\5, 135. 
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Figure 19. Reflection coeffilcients {amplitude and 
phase) of the lower ionosphere for low frequencies 
during blackout conditions illustrating the fre- 
quency dependence of the classical magneto-ionic 
theory. 

\Tee\ and | Tmm\ arc also illustrated for comparison with auiescent 
conditions (see fig. 5-9 for complete qUiCSCcnt reflection coefficients). 
0a=225, 315. 
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Figure 20. Reflection coefficients {arnplitude and 
phase) of the lower ionosphere for low frequencies 
ilhistrating the action of the Lorentzian force {earth's 
magnetic field) in the classical magneto-ionic theory 
during blackout conditions. 

1 Tee\ also shown for comparison with auiescent conditions. (See 
fig. 10-14 for complete quiescent reflection coefTicients.) /=10 kc/s. 



6. Conclusions 

The reflection coeflScients of the lower ionosphere can be more rigorously determined from 
geophysical data on the lower ionosphere with the aid of techniques developed in this paper. 
The classical magneto-ionic theory can be employed, but is subject to restrictions as to the 
use of the average collision frequency, g^v. It is therefore concluded that the same problem 
should be re-investigated with the aid of the monoenergetic electron collision frequency and the 
formulation presented in this paper for the determination of the complex parameter, g] com- 
parison can then be made with the classical theory. 
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FiGUEE 2i. Reflection coefficients {amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the action of the Lor entzian force {earth's 
magnetic field) in the classical magneto-ionic theory 
during blackout conditions. 

I Tee\ also shown for comparison with quiescent conditions. (See 
fig. 10-14 for complete quiescent reflection coefficients.) /=20 kc/s. 
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Figure 22. Reflection coefficients {amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the action of the Lorentzian force {earths 
magnetic field) in the classical magneto-ionic theory 
during blackout conditions. 

\ Tee\ also shown for comparison with quiescent conditions. (See 
fiR. 10-14 for complete quiescent reflection coefficients.) f=50 kc/s. 
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Figure 23. Reflection coefficients {amplitude and 
phase) of the lower ionosphere for low frequencies 
illustrating the action of the Lorentzian force {earths 
magnetic flM) in the classical magneto-ionic theory 
during blackoid conditions. 

I Tee\ also shown for comparison with Quiescent conditions. (See 
fig. 10-14 for complete quiescent reflection coefficients.) ;=100 kc/s. 
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Figure 24. Reflection coefficients {amplitude any 
phase) of the lower ionsophere for low frequencies 
illustrating the action of the Lorentzian force {earths 
magnetic field) in the classical magneto-ionic theory 
during blackout conditions. 

I Tee\ also shown for comparison with quiescent conditions. (See 
fig. 10-14 for complete quiescent reflection coefficients.) /=200 kc/s. 
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The results of this investigation indicate that the total blackout phenomena which 
characterize HF do not extend through LF and VLF for the model ionosphere investigated 
with the classical theory. Indeed, for the situations computed here the field was found to be 
enhanced at VLF. The emplo}Tiient of these models does not alter previous conclusions on 
the action of the Lorentz force on the electron in the ionosphere, i.e., the propagation mto the 
cast is enhanced with respect to the propagation into the west, except to a lesser degree at 
the higher frequencies, which also exhibit increased attenuation. 
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